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RADIOFREQUENCY SURFACE DETECTION COIL 

Technical Field 

The present invention relates to a coil for detecting a magnetic resonance signal induced in a sample and, 
more particularly, to a noise-immune coil having an optimized magnetic field strength in a specified inspection 
region outside the coil. 

Background Art 

Magnetic resonance is useful to detect the presence of a specific substance in a sample. For example, 
generally, radio frequency (RF) radiation at a particular frequency will induce a magnetic resonance signal in a 
specific substance, but not in other substances. Therefore, the induced magnetic resonance signal can be detected 
to thereby indicate the presence of the specific substance. 

It is common to detect a magnetic resonance signal by placing a sample to be measured in or near a tuned, 
electronically resonant tank circuit. Then, the response of the tank circuit to the electromotive force produced by 
nuclear or electronic spins in the sample is measured. With Nuclear Magnetic Resonance (NMR), Nuclear 
Quadrupole Resonance (NQR), and low-frequency Electron Paramagnetic Resonance (EPR) the sample is placed 
in or near an inductor, commonly referred to as a coil, that detects AC magnetic fields. The inductance of the coil 
is tuned with a parallel and/or series capacitance to make the circuit electrically resonant at the measurement 
frequency. One or more additional reactive impedances (inductors or capacitors) are typically added to adjust the 
resistive impedance at resonance to a particular value which optimizes the detection sensitivity. 

NQR detection systems for the detection of explosives and narcotics, and various NQR concepts, are 
disclosed, for example, in U.S. Pat. No. 5,233,300, "Detection Of Explosive And Narcotics By Low Power Large 
Sample Volume Nuclear Quadrupole Resonance (NQR)"; U.S. Pat. No. 5,365,171, "Removing The Effects Of 
Acoustic Ringing And Reducing Temperature Effects In The Detection Of Explosives By NQR"; U.S. Pat. No. 
5,206,592, "Detection Of Explosives By Nuclear Quadrupole Resonance"; U.S. Pat. No. 5,608,321, "Method And 
Apparatus For Detecting Target Species Having Quadrupolar Nuclei By Stochastic Nuclear Quadrupole Resonance"; 
and U.S. Pat. No. 6,054,856, "Coil Which Is Immune To Environmental Noise", all of which are incorporated herein 
by reference. 

In most NQR, NMR and EPR applications, a common coil would typically be used as both a receiver coil 
and a transmitter coil, though this is not essential. In MRI separate coils would typically be used. The asymmetric 
gradiometer coils can be used for either or both of the transmitter and receiver coils. 

Figure 1 is a diagram illustrating an example of a conventional magnetic resonance apparatus. Referring 
now to Figure 1, a transmitter 20 and a receiver 22 are connected to a probe 24 through a transmit/receive (T/R) 
switch 26. Probe 24 includes a coil 28, forming part of a resonant, tuned tank circuit 27 with various other inductors 
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L and capacitors C as tuning elements. To detect the presence of a target substance, T/R switch 26 connects 
transmitter 20 to probe 24 while disconnecting receiver 22 from probe 24. Then, transmitter 20 generates a pulse 
and supplies the pulse to probe 24. As an example, in NQR, the pulse is formed from an RF signal having a 
frequency corresponding to the resonance signal of the target substance which is intended to be detected. Probe 24 
receives the pulse, which causes coil 28 to store (RF) energy. 

If a sample (not illustrated) is appropriately placed near or inside coil 28, the stored RF energy will cause 
a corresponding RF magnetic field to irradiate the sample. If the sample includes the target substance, the RF 
magnetic field will induce a magnetic resonance signal in the target substance. For example, if the apparatus 
operates under the principles of NMR (which includes Magnetic Resonance Imaging, MRI), then an appropriate 
NMR resonance signal will be induced. If the apparatus operates under the principles of NQR, then an appropriate 

NQR resonance signal will be induced. 

After the sample is irradiated with the RF magnetic field, T/R switch 26 connects receiver 22 to probe 24 
while disconnecting transmitter 20 from probe 24. Coil 28 then detects the resonance induced in the target 
substance, and probe 24 produces a corresponding output signal. The output signal of probe 24 is received and 
analyzed by receiver 22, to confirm the presence and/or measure the quantity of the target substance in the sample. 

Figure 1 is only one example of a magnetic resonance apparatus. For example,Figure 1 illustrates T/R 
switch 26 to connect transmitter 20 and receiver 22 to the same probe 24. However, instead, a transmitter and 
receiver can each have a separate, dedicated probe together with a switch or gate for protecting the receiver while 
the transmitter is ON. Note that such an apparatus may also be used to measure material specific RF signals which 
arise from electromagnetic properties other than the magnetic resonance phenomenon. 

Figure 2 is a diagram illustrating a simple, conventional coil which can be used in a probe. Referring now 
to Figure 2, a coil 29 typically forms a loop. Typically, a tuning capacitance C and a matching capacitance C are 
also provided. The coil shown has 1 turn, though additional turns are often used. 

In magnetic resonance the signal-to-noise ratio (SNR) is determined, in part, by the noise contributions and 
the quality factor (Q) of the receiver coil. It is well-known that random thermal noise contributions typically arise 
from Johnson noise in the RF inspection coil and the first amplifier in the receiver. In the case where the probe 
cannot be electrically shielded, a further noise contribution arises from extraneous environmental noise. 

It is also well-known in magnetic resonance that the Q of a receiver coil is determined not only by resistive 
loss in the windings of coil itself but also by loss in nearby electrically conducting materials that can dissipate energy 
from currents in the coil. As the SNR typically varies as Q*\ such electrical loss in the surroundings leads to a 
reduction in SNR. For example, in MRI, the main source of electrical loss can come from the patient, and not the 
receiver coil windings, as the water in the body has an electrical conductivity comparable to sea water. In NQR land 
mine detection it is found that some soils also present significant electrical loading to the receiver coil, leading to 

a decreased coil Q and decreased SNR. 

In many applications in MRI and also in land mine detection, a surface coil is used for the inspection. The 
larger the surface coil, the more the receiver Q is decreased by the presence of electrically conducting materials in 
the surroundings. Indeed, in MRI, a system designer conventionally chooses as the receiver coil the smallest surface 
coil that will "cover" the region of interest. 
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For a conventional simple circular coil (Figure 2) of radius R the RF magnetic field falls off with distance 
so the coil is very inefficient at receiving NMR or NQR signals much beyond a distance R. Hence, for applications 
such as MRI or land mine detection, the distance to the inspection or interrogation region of interest determines the 
minimum size of the coil which is useful. For such a circular coil on the surface of a large conducting volume, we 
5 find that approximately 65% of the electrical loss from this volume arises from regions that are deeper than a 

distance R below the surface. Hence, most of the loss comes from a region beyond the actual region that can be 
imaged (in MRI) or inspected (in NQR landmine detection). 

Although an RF coil used as a detector in magnetic resonance is not specifically designed to detect or 
receive radio signals, the detection coil does tend to pick up such unwanted interference. This radio interference from 
10 radio stations or other RF noise sources in the relevant frequency range can overwhelm the magnetic resonance 

signals of interest. One solution is to employ an external RF shield surrounding the coil and sample, however, this 
is impractical in many applications and can unduly increase the cost and size of the system. For the simple circular 
coil, the susceptibility for receiving interference from distance RF sources increases in proportion to the square of 
the coil radius, R. 

1 5 Accordingly, it is desired to reduce the environmental noise pickup and reduce the electrical loss due to 

the proximity of a conducting medium while at the same time keeping adequate sensitivity to the desired signal. 
Approaches for achieving this may employ detector (or transmitter) coils that are designed to (i) strongly reject 
electromagnetic environmental noise from more distant sources, (ii) to couple strongly to the desired signals that 
arise from the desired inspection region very near the coil, and (iii) to couple more weakly to other nearby 

20 electrically conducting materials. 

In most NQR, NMR and low-frequency EPR applications, a common coil would typically be used as both 
a receiver coil and a transmitter coil, though this is not essential. In MRI separate coils would typically be used. 
By the principle of reciprocity, all coil designs can be evaluated by considering them for use as a transmitter even 
if the coil is ultimately intended to be used as a receiver. That is, if when used as a transmitter a coil produces a large 

25 magnetic field at a specified location, that same coil can be expected to be a good receiver for magnetic resonance 

signals arising from material at that location. 

One approach to reduce interference is 'balancing', in which the interference enters but is arranged to cancel 
itself based on some property of the detection coil and the external interference. A gradiometer is one example of 
this approach, as it can be designed to respond to any order of the spatial derivative of the electric or magnetic field. 

30 A conventional linear or first order gradiometer is primarily sensitive to the first spatial derivative of the 

field, and is correspondingly insensitive to field components that do not vary in space. Since the spatial variations 
of an electromagnetic wave are characterized by its wavelength, this provides a means to cancel fields that have 
wavelengths much larger than the characteristic size of the gradiometer, thereby reducing noise pickup from distant 
sources. (For example, the free space wavelength of a 1 MHz signal is 300 meters whereas the coils will generally 

35 be less than 1 m across.). A simple conventional (linear) magnetic field gradiometer is formed of two loops which 

are spatially removed from each other and have currents flowing in opposite directions. 

Figure 3 is a diagram illustrating a conventional magnetic field gradiometer. Referring now to Figure 3, 
a conventional magnetic field gradiometer has a conductor 30 which forms two loops 32 and 34 of identical radii 
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but wound in the opposite sense. Each of these loops may have more than 1 turn of wire. The direction of the 
current changes in coils 32 and 34, so that magnetic fields generated by coils 32 and 34 are opposite each other. 

For an application such as buried land mine detection, the inspection volume is located in one half space 
(underground) and the detection apparatus is limited to the other half space (above ground). The coil is used as part 
5 of a tuned circuit to detect the NQR signal which for landmine detection is a narrowband (~ 1 kHz) signal in the range 

of 0.5-5 MHz. These NQR signals are very weak and RF interference is a major problem for such field work. 
While a conventional magnetic field gradiometer will reduce environmental magnetic field noise pick-up, the amount 
of reduction is often insufficient for many applications. In addition, there is also a reduction of the desired signal, 
compared to the simple circular coil, which is larger than desired. Here the use of an external RF shield, which must 
10 surround both the coil and the land mine, is not practical. It is therefore desirable to provide an improved 

gradiometer-type coil design having improved immunity to interference without extra shielding and at the same time 
providing an efficient detection region outside (i.e. near, but away from) the coil. 

Disclosure of the Invention 

15 According to the invention, a magnetic field gradiometer detector for detecting a material of interest 

positioned in a detecting region outside the detector includes a transmitter for generating an output signal at a 
selected frequency, a receiver for detecting the signal, a probe, and a switch coupling the receiver and transmitter 
with the probe that alternately connects and disconnects the receiver and transmitter to the probe, switching between 
transmitting and receiving. The probe includes tuning elements and a gradiometer coil array. The gradiometer coil 

20 array includes a first surface coil and a second surface coil wound in an opposite sense, the probe having a first side 

and an opposite second side, with the first and second surface coils configured asymmetrically such that the probe 
projects a magnetic field in the outside detecting region adjacent to the first side while being self-shielded on the 
second side. Accordingly, the material of interest is detected with increased sensitivity and accuracy in the detecting 
region with the detector shielded from outside R sources. 

25 Also according to the invention, a method for analyzing a sample by magnetic resonance includes the steps 

of (a) generating a train of radio frequency pulses having a predetermined frequency, (b) transmitting the train of 
radio frequency pulses to the magnetic gradiometer coil array, (c) irradiating the sample in response to the train of 
radio frequency pulses transmitted to the coil at step (b), (d) detecting a signal in response to irradiating the specimen 
at step (c), and (e) receiving the signal detected at step (d). 

30 The invention provides improved immunity to interference from outside RF sources, reduced losses due 

to nearby electrically conducting materials, and an improved detection region outside the coil, thereby providing 
increased sensitivity for detecting small quantities of target materials. These properties are advantageous in many 
applications, such as in buried land mine detection and humanitarian de-mining, security inspections, contraband 
detection, and in geophysical surveying activities such as well logging. 

35 Additional advantages of the invention will be set forth in part in the description which follows, and, in part, 

will be obvious from the description, or may be learned by practice of the invention. 
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Brief Description of the Drawings 

These and other objects and advantages of the invention will become apparent and more readily appreciated 
from the following description of the preferred embodiments, taken in conjunction with the accompanying drawings 
of which: 

5 FIG. 1 (prior art) is a diagram illustrating an example of a conventional magnetic resonance apparatus. 

FIG. 2 (prior art) is a diagram illustrating a simple, conventional surface coil. 
FIG. 3 (prior art) is a diagram illustrating a conventional magnetic field gradiometer. 
FIG. 4 is a diagram illustrating a basic coil geometry for use in deriving optimal coil parameters according 
to the invention. 

10 FIG. 5 is a graph showing the magnitudes of the ideal current densities computed using parameter /and 

a single-layer surface coil for two inspection vokimes using the corresponding eigenfunctions. 

FIG. 6 is a graph showing the magnitudes of the ideal current densities computed using parameter g and 
a single-layer surface coil for two inspection volumes using the corresponding eigenfunctions. 

FIG. 7 is a diagram illustrating ideal current densities that maximize /for one inspection volume for a two- 
1 5 layer gradiometer surface coil according to an additional embodiment of the present invention. The solid line is the 

current in the layer that is near the inspection region, the dashed line in the layer which is away from the inspection 
region (the "shield layer"). 

FIG. 8(a) is a schematic diagram illustrating a method to approximate the ideal current densities in one 
layer using unequally spaced circular windings. 
20 FIG. 8(b) is a schematic diagram illustrating the connections for a two-layer concentric loop coil, each layer 

being wound in a manner similar to that of FIG. 8(a), though with different numbers of windings and different 
spacings, and showing the coil interrogation region according to an embodiment of the present invention. 

FIG. 9(a) is a diagram illustrating a cross section of a two-layer concentric loop coil of FIG. 8(b), though 
with more windings in each layer, according to an embodiment of the present invention. 
25 FIG. 9(b) is a diagram illustrating a cross section of a two-layer coil using wide flat wires according to an 

embodiment of the present invention. 

FIG. 10 is a diagram illustrating the results of the finite element computations for four self-shielded 
gradiometer configurations according to an embodiment of the present invention. Figures labeled c and d 
correspond to the coils shown in figures 9A and 9B respectively. 
30 FIG. 11 is a diagram showing a comparison of the magnetic field profiles along the coil axis (r = 0, z) for 

the computed field and the corresponding measured values for one of the prototype coils according to an 
embodiment of the present invention. 

FIG. 12 is a diagram showing the computed field (a) which would be produced by one unit of current using 
only the bottom coil and the field (b) from the two-layer shielded coil according to an embodiment of the present 
35 invention. 

FIG. 13 is a diagram illustrating computed values for the on-axis magnetic fields for the two cases of FIG. 
12 according to an embodiment of the present invention. 
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FIG. 14 is a graph showing the expected variation in SNR for a particular inspection region and for 
different interference levels and coil layer spacings according to an embodiment of the present invention. 

Best Mode for Carrying Out the Invention 

5 Reference will now be made in detail to the present preferred embodiments of the present invention, 

examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to like 
elements throughout. 

Definition(s): As used herein, the term "asymmetric" when referring to the plural coils of a detector or coil 
array means that the coils are non-identical, that is, at least one coil of the plurality of coils has a different geometry, 
10 different number of turns or windings, or differs in some other physical or operating parameter from the other coil 

or coils of that embodiment. 

Theory used to derive coil designs 

For our purposes it is sufficient to assume that the RF frequency is low enough that displacement current 

15 terms can be neglected -that is, we are in the long wavelength limit. For NQR frequencies of the order of 1 MHz 

and coil dimensions of 1 m or less, this is a reasonable assumption. In addition, for this derivation we assume all 
regions outside the coils are free from magnetic, dielectric, and electrically conducting media and we specialize to 
the cylindrically symmetric problem with planar surface coils. We also will assume that the coil is designed and 
constructed so that it is sensitive only to magnetic fields and not directly to electric fields. 

20 We first examine a single-layer coil and develop a computational methodology for optimizing the current 

distribution using a variational approach and/or a diagonalization of the corresponding eigenvalue problem. We then 
extend to the case of a gradiometer using two coaxial coils and demonstrate both theoretically and experimentally 
that such a coil can be constructed to be self-shielded in one half space but still produce a usable field in the other 
half space. We then consider the trade-off between SNR and immunity to RF interference for such a gradiometer 

25 when used for NQR detection. 

Single-layer coils 

We consider the basic coil geometry shown in Figure 4. The inspection volume of interest is taken to be 
a cylinder of radius a in the half-space z > 0 . The coil is in the plane z = 0. The ends of the cylinder are at z = b 
30 and z = c , b < c, as shown. The axis of the cylinder defines the z-axis of our cylindrical coordinate system. For the 

computations included here, units are chosen so that a - 1. Our main goal is to maximize the field in the inspection 
region which is produced by the coil when used as a transmitter, while at the same time minimizing the field 
produced elsewhere. 

In the region away from the coil, z *0, it is straightforward to show that the general solution to Maxwell's 
35 equations in cylindrical coordinates, neglecting the displacement current term and restricted to finite solutions with 

cylindrical symmetry, can be written 
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oo 



B z = \aA(a)J 0 (ar)e' a[zl da 

o 

B r = - j aA(a)J^(ar)e' aM da = j aA(a)J x (ar)e' aM da 

0 o 

B 0 = O 



where J 0 and J } are Bessel functions and the "expansion coefficients," A(a), are at this point arbitrary. The 
boundary condition at z = 0 provides the relation between the expansion coefficients and the surface current density, 
K c (r), 



K r (r) oc zx B(r,z = 0) = 5 r (r,z = 0) 3 = £ j 054(a) J x (at)da (2) 

o 

where the proportionality constant depends on the units used. 

There are several useful quantities we will compute starting with equation 1 . For example, the total energy 
in the magnetic field (in convenient units) is 



E tot = \d 3 x\B\ 2 = 47i\~ da A 1 (a) . 



(3) 



all space 



For coils with uniform resistivity, the dissipation in our coil, W, is 



Woe ^IjtrdrK] = 2n\daaA 2 (a) (4) 



0 0 



10 and the magnetic energy within the cylindrical inspection volume is 



(2 C oo oo 



E v = 1n\ rdr\dz \B\ 2 = 1n\ da\ da! A(a')G(a' ,a)A(a) , (5) 



Ob 0 0 



where 



7 



Inventor(s): Allen N. Garroway et al. 



Patent Application 

Navy Case No. 84,294 



10 



G(a',a)=a %^ {e^ b - e^ a ' c )(j,{aa)J 0 (a'a)- J 0 (aa)J x {a'aj) ( 6 > 

and the proportionality constants for E m and E v are the same. 

It is useful to consider a discretized version of the above equations. That is, one considers N discrete values 
of a with spacing Act. Then the expansion coefficients A(0C) become a vector, A, of length N, and G(a',a) becomes 
an NxN symmetric real tensor, G. Since N can be arbitrarily large and Ace arbitrarily small, for practical purposes 
this discretization can be considered a change only in notation. Later, we will use specific values for N and Aa, in 
which case one should verify that A and G adequately represent the continuous functions in the regions where those 
functions are non-zero. Typically this verification is done by making a significant change to the sampling grid, N 
and/or Aa, and verifying that the solution does not change significantly. With this new notation, equations 3-5 
become 

E tot = A T -A 

W = A ,T .A' (7) 
£ V = A T -G-A , 



where A' = -yf(X A . 



The problem now is to determine an appropriate choice for A, from which the current densities can be 
computed using equation 2 . We now consider several different criteria for choosing A. 

One criterion to optimize is the "filling factor." While there is no precise definition of filling factor for the 
15 general surface coil case, we still wish to characterize the fraction of the total magnetic field which is in the 

inspection volume. Ideally this would be 100% so that all of the RF field produced would be useful. For our 
infinite region, we can define a filling factor,/, as the ratio of the energy in the magnetic field within the inspection 
volume to the total magnetic field energy everywhere. That is, 

F A T • G A 

f = ^= A ,p A . (8) 

1 E tot A T • A 



20 Since G is a symmetric real tensor, it can be diagonalized. The maximum value of/can thus be determined by 

finding the largest eigenvalue of G. The appropriate choice of A will then be the eigenvector corresponding to that 
largest eigenvalue. 
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There are two approaches we have looked at to determine A which produces the maximum value of /- a 
variational approached based on a numerical implementation of the continuous integral equations, and a numerical 
diagonalization of the discrete tensor, G, using standard numerical subroutines. Other methods can be used as well. 

For the variational approach, a trial function with one or more variable parameters is used for A. The 
parameter(s) is (are) varied numerically to search for a maximum value of/. The trail function which gives the 
maximum value is then an approximation to the desired result. The method is time consuming and it may be 
necessary to test several distinctly different trial functions to achieve accurate results, but the variational approach 
is useful since it will provide a measure of the sensitivity to small changes in the function A(ct), and hence the 
sensitivity to changes in the current distribution. In addition, there is considerable flexibility for choosing trial 
functions that match known physical situations which can be used for comparison. 

Since G(a\a) falls off rapidly when |<x'-<x| > l/a , a numerical diagonalization of G proceeds rather 
quickly. However, if N is made too large, numerical round-off problems may arise. We have used N as large as 200 
with good results. The numerical diagonalization can be done on a desktop personal computer in less than 1 second. 
Figure 5 shows the magnitudes of the current densities computed from the corresponding eigenfunctions for two 
different interrogation regions. The resulting values off compare well with the largest values found using the 
variational technique and hence either technique can be used in practice. 

Other parameters can be maximized in the same way. For example, to maximize the magnetic field energy 
in the inspection volume for a given power dissipation in the coil, the "excitation efficiency," we maximize 

A' T - G'.A' (9) 
s W A • A' 



where 

G'(a\a) = G(a\a)/yfaa 7 (10) 

which has the same mathematical form as was used for the filling factor. The procedure is then to diagonalize G' 
and find the eigenvector A ' that has the largest eigenvalue. Then one computes A from A ' and A is used to find the 
desired current densities as before. The current densities that maximize g (for the two inspection volumes used 
previously) are shown in Figure 6. These current densities are very similar to those that maximize the filling factor, 
except the tail of the current density extends to somewhat larger r. The fact that these are similar means that a coil 
which approximates either of the two continuous distributions shown will be close to optimum for both of these, and 
other similar, criteria. For example, if one uses the results which optimize g,/is reduced by only about 10 to 15% 
and vice versa. 

Above we have used as examples the two parameters / and g for the optimization. Of course other 
parameters, or modifications to these parameters, may be more appropriate for specific applications and are easily 
derived by those practiced in the art. 

We do point out that not all parameters yield useful results. For example, the parameter q, defined by 
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E 

is related to the quality factor of the coil. We find that when q is maximized, the resulting current density is nonzero 
only at infinite radius and that for finite currents the field produced is infinitesimally small. A similar aberrant 
solution results when one attempts to maximize the homogeneity of the magnetic field. In the limit, the field can be 
5 made exactly uniform (though vanishingly small) over the inspection volume by pushing the current density out to 

infinitely large r. 

■t 

Two-layer coils - the self-shielded axial gradiometer 

The geometry we consider here is identical to that for the single-layer coil except a second layer that may 
10 contain current is added at z = -d . Using subscripts 1 and 2 to represent the layer at z = 0 and z = -d respectively, 

we write the total field as the superposition of the fields from the currents in the individual layers. That is, for 
example, 

oo 

B z = j daaJ 0 (cxr)(A l (a)e' a ^+ A 2 (a)e~ alz+<il ) . (12) 

o 

We can hope to achieve good noise immunity at low frequencies if the field produced by this coil is 
identically zero for z < -d, sometimes referred to as "supershielding ." By reciprocity, the coil will not receive 
15 signals from interfering sources in that region. At RF frequencies and for distant sources such a solution is not 

easily achieved in practice, however the radiated energy - the "quadrupole radiation" - will be smaller than that of 
a comparably sized single-layer coil by a factor roughly equal to the ratio of the coil spacing to the wavelength. In 
our case the wavelengths range from 30 to 300 m while a representative spacing is 10 to 100 cm, and so this 
reduction can be quite substantial. This reduction is achieved with a configuration for which at dc (zero frequency, 
20 infinite wavelength) the supershielding condition is satisfied. 

Supershielding (at zero frequency) will be achieved if 

A x (a)= - A 2 (a)e ad (13) 

and with this condition, the field in the region z > 0 (where the inspection volume is located) has the same 
mathematical form found for the single-layer coil, 

oo 

B z (z>0)= \daa A eff (a) J 0 (ar)e~ az (14) 

o 

10 



(11) 
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where 

-2ad 



A e# = A,(l+e^) (15) 



and there is a similar expression for B r . The total magnetic energy is 



E tot = \ 27rrdr\dz\B\ 2 = An\ da A 2 (a)(l- {2ad + l)e" 2arf ) 



0 -~ 0 

00 

2 



= An\ da(A(a)) 



where 



(16) 



A(a ) = A, (a )y/l-(2ad+ l)e' 2ad (17) 



and so the filling factor can then be written as before 

AL. G a~ a t g a 

r _ g ff w_ _ _ 2l (18) 

7 A -A A -A 

where 

Vl" (2ad + l)e' 2ad Vl- (2tf'rf + l)e" 2arf 



which is once again an eigenvalue problem. The largest eigenvalue of G and corresponding eigenvector, A , 

are found. Then by working backwards, A x and A 2 and the corresponding current densities in the two layers are 
determined. Figure 7 shows the results where / has been optimized for one particular inspection volume and 
spacing between the layers. The current in layer 2 is reduced compared to layer 1, and is pushed out to a larger 
10 radius - that is, the gradiometer produced is asymmetric. The negative contribution to the magnetic field within 
the interrogation region from layer 2 is significantly reduced compared to a symmetric gradiometer (A, = -A 2 ) with 
the same spacing. The net magnetic dipole moment, however, remains zero, at least within the accuracy of a 
numerical calculation. It is interesting to note that the filling factor,/, as defined above, for this two-layer coil is 
actually 30% larger than that of the corresponding optimized single layer coil. 
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A similar procedure would be used to derive ideal current densities for the parameter g, or other similar 
parameters. 

Practical Implementation 

5 For the design of real coils that produce a magnetic field which approximates the field found using the 

continuous theory above, there are several practical considerations. One obvious limitation is that the coil should 
be finite in size. For use under field conditions the coil assembly must ideally be small and lightweight, especially 
for a hand-carried device. For the cases we have considered, with the exceptions noted above, the current densities 
fall rapidly to zero with increasing radius and so in principle, size should not be a major problem. 

10 In addition, while it is possible to produce a finite continuous sheet with variable resistivity in order to 

obtain a continuously variable current density, that is not an appropriate procedure here. The main goal here is to 
achieve a sensitive RF receiver coil. Adding resistance to the coil to control the current will, however, add noise 
which will subvert our attempts at making a sensitive receiver coil. That is, we desire approximate solutions that 
do not also unduly compromise the coil quality factor, Q. 

15 To approximate the continuous current density without severely compromising Q we use an array 36 of 

concentric circular loops 38 and 40 wired in series but with unequal spacings, as shown in Figures 8(a), 8(b, 9(a), 
and 9(b). Current is brought into and out of each loop 38 and 40 in the radial direction. The leads 42 and 44 are 
closely spaced (the spacing is exaggerated in the figure for clarity) and run in the radial direction and can be thought 
of as a short length of transmission line and will contribute an insignificant amount to the magnetic field in the 

20 inspection region and at larger distances. The spacing between the individual circular loops in each layer 38 and 

40 is proportional to the inverse of the computed current density. Accordingly, coils 38 and 40 are stacked and wired 
in series, but with current flowing in the opposite sense. It is also possible to use other means than direct coupling 
through wires and a series connection between the coils to create the desired currents, such as via various capacitive 
or inductive couplings, as is well-known by those practiced in related art. It should also be noted that coils of other 

25 shapes (polygonal, etc.) can also be used to approximate the desired current distributions, though possibly less well, 

and may be more convenient in some cases. 

Here we present the results for two prototype two-layer coils for comparison to theory. One of these coil 
arrays, coil array 46 shown in Figure 9(a), uses a relatively large number of concentric loops to provide a reasonably 
accurate approximation to the continuous results. The other, a more practical design, is a fairly crude approximation 

30 using a small number of wide (ca. 2.5 cm) flat wires cut from copper sheet, shown in Figure 9(b). 

Magnetic field profiles were calculated taking advantage of the axial symmetry and using a two-dimensional 
electromagnetic finite element simulation software package (Ansoft Corporation Maxwell 2D version 9) . These 
simulations do not include the wires leading into or out of the loops. Typically 20K or more finite elements were 
used. The simulations were performed at both dc and RF frequencies (~3 MHz) using copper conductors, 

35 appropriate for landmine detection: little difference could be seen. 

Figure 10 shows results of the finite element computations for four self-shielded gradiometer configurations 
starting with the continuous current distribution (see Figure 7), an approximation to the continuous distribution 
using a large number of discrete wires carried out to large radius, the same discrete wires but truncated at a smaller 
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radius, and a coil constructed using a relatively small number of discrete, but wide conductors, respectively. The 
latter two represent two prototype coils used for measurements. It can be see that the magnetic field profiles (away 
from the discrete wires) are very similar for all of these. 

Figure 11 shows a comparison of the magnetic field profiles along the axis (r = 0, z) for the computed field 
5 and the corresponding measured values for one of the prototype coils. The measurements were made at RF 

frequencies (3 MHz) using a 5 cm diameter, 3-turn pick-up coil and have been normalized to the value at z = a. 
For the numerical results the flux through a 5 cm diameter pickup coil was computed assuming the pick-up coil had 
no depth along the z-axis. The theory and experiment agree reasonably well. Experimentally, the measured signal 
does not decrease as rapidly on the inspection side of the coil as the computations suggest. 
10 The addition of the second layer to provide shielding will have an undesirable impact on the magnitude of 

the field in the inspection region. Figures 12 and 13 illustrate this effect for one of our two-layer coils. In Figure 
12, (a) shows the computed field which would be produced by one unit of current if only the bottom coil is used (this 
field profile is essentially identical to what is obtained for the optimized single-layer coil), while (b) shows the field 
from the two-layer shielded coil (with d = a) when the same current is used in the bottom coil. Computed values 
1 5 for the on-axis magnetic fields for these two cases are shown in Figure 13 where the self-shielding on the side away 

from the sample (z < 0) and only a small reduction in the field in the inspection region ( i < z/a < 2 ) can be seen. 
It is also pointed out that the reduction in the field due to the second layer is less than that obtained for an equivalent 
conventional or "anti-symmetric" gradiometer with the same layer spacing, e.g. obtained using Aj(ct) = -A 2 {Ct). The 
improvement over the anti-symmetric case depends on the layer spacing; it is about 15% near the center of the 
20 inspection region for the case shown in Figure 10. 

To achieve less reduction of the field in the inspection region, the two layers can be moved farther apart, 
though the diameter of the shield layer will also need to increase, with a corresponding decrease in the net current, 
to maintain the same degree of shielding. At the same time, however, the immunity to distant RF sources will be 
compromised as the spacing becomes a larger fraction of the wavelength, as discussed above. To determine an 
25 optimum layer spacing will depend on the size of the interfering RF signal(s). With no RF interference, an infinite 

layer spacing (i.e. removal of the shielding coil) yields the largest SNR. 

To determine the optimum gradiometer spacing in detail requires knowledge of the changes in the size of 
the RF quadrupole radiation produced by the coil when used as a transmitter, the direction(s) to the source(s), and 
the expected size of the interference compared to thermal noise. That is, the optimum spacing will depend upon the 
30 ultimate application. Such a complete calculation is beyond the scope of this work. However, an estimate of the 

appropriate spacing can be made for relative changes one might expect. 

The gradiometer is sensitive to the spatial gradients of the magnetic field measured over the size scale of 
the coil. Hence, to a reasonable approximation one can infer that for our gradiometer coils the sensitivity to radiation 
from distant sources will be proportional to the layer spacing, d. For NQR detection in the absence of interference, 
35 the SNR for the near-field NQR signal will depend on the product of two terms: the square root of the coil quality 

factor, Q, and the size of the magnetic field produced by the coil in the inspection region for a unit current excitation. 
We then proceed as follows. The optimum continuous current densities are computed as described above. In the 
example shown here the parameter/is being optimized. Using those continuous current densities, we compute the 
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magnetic field within the inspection region as a function of spacing, B(d). In our case we did this latter step using 
the finite element software discussed above after normalizing the currents so that the total circulating current in the 
lower coil was unity. The field at the center of the inspection region was used. The changes in the coil quality factor 
with separation, q(d), are estimated using the results of the optimization and equation 13. Then the dependence of 
the SNR on the gradiometer spacing and in the absence of interference, SNR 0 (d) is 



w i^um (20) 



where d 0 is a spacing chosen as a reference point for comparison. Equation 20 then represents the size of the signal 

compared to thermal noise. 

We assume the interfering signal has a coherence time long compared to Ql V, where V is the operating 
frequency. This type of narrow band signal is typical of interference observed from broadcast stations. In that case 
10 the observed interfering signal also is proportional to the square root of the quality factor. Defining the parameter 

a as the magnitude of the interfering signal for a reference gradiometer spacing d 0 expressed in units where the 
thermal noise is 1, and using the approximation that the sensitivity to interference is linear in d, we then find SNRj 
, the ratio of the desired signal to the sum of the thermal noise and the interfering signal, will depend on d according 
to 

SNR 0 (d) 

SNR i ~ ~, — TT^T • (21) 

' 1 + ad I d 0 

15 For the two-layer gradiometer, Figure 14 shows the results for various values of df, expressed in dB, one particular 

inspection region, and d 0 ~6a . 

Hence, for the inspection region used here and when interference is large (e.g. a= 30 dB), a layer spacing 
of approximately a, the inspection region radius, should be used as seen in Figure 14: for that spacing the NQR 
SNR is reduced to about 30% of what one could achieve in the absence of interference. To achieve some limited 
20 degree of interference reduction without significantly reducing the NQR sensitivity, one would choose a spacing of 

about 3a or more. In that case some other means might be employed to separate the remaining interference from 
the desired signal. 

The invention is relatively easy to implement in practice using approximations, as described above, to 
results computed using either a variational approach or through standard numerical matrix diagonalization 
25 techniques. The technique is applied in particular to develop a practical self-shielded asymmetric gradiometer that 

is relatively immune to RF interference from distant sources with a minimal compromise to its ability to sense RF 
from nearby sources, properties necessary for NQR detection of buried land mines. 

According to the above embodiments of the present invention, the coil is virtually immune from common 
environmental noise sources, both electrical and magnetic, and hence eliminates the necessity of providing external 
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RF shielding. This will have utility for NQR detection of explosives and perhaps obviate the need for external RF 
shielding in magnetic resonance imaging (MRI) applications. In particular, embodiments of the present invention 
can be applied to the use of NQR to the detection of nitrogenous or chlorine-containing explosives or narcotics 
carried in luggage, mail, small cargo, on a person, or buried in the ground. More generally, the above embodiments 
of the present invention can be used for improved detection of any material which has a magnetic resonance response 
that is unique compared to its surroundings. 

Moreover, the above embodiments of the present invention are applicable to many different areas of RF 
signal detection. For example, the above embodiments of the present invention are applicable to NMR, NQR, MRI, 
broad line NMR studies such as stray field imaging of materials and possibly to pulsed Electron Paramagnetic 
Resonance (EPR) spectrometers at lower frequencies. Therefore, for example, the embodiments of the present 
invention are applicable to the detection of NMR resonance signals, NQR resonance signals and EPR resonance 
signals. 

In most NQR, NMR and EPR applications, a common coil would typically be used as both a receiver coil 
and a transmitter coil, though this is not essential. In MRI separate coils would typically be used. The asymmetric 
gradiometer coils can be used for either or both of the transmitter and receiver coils. 

There are many different types of probe configurations which can use a coil according to the above 
embodiments of the present invention. Thus, the present invention is not intended to be limited to any specific probe 
configuration. 

Although a few preferred embodiments of the present invention have been shown and described, it would 
be appreciated by those skilled in the art that changes may be made in these embodiments without departing from 
the principles and spirit of the invention, the scope of which is defined in the claims and their equivalents. 



15 



